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1 .  Introduction 


Consideration  of  long  wavelength  infrared  radiation  emitted  by  the 
terrestrial  atmosphere  is  important  in  application  of  remote  sounding  to  the 
atmosphere,  treatment  of  the  heat  budget  of  the  planet,  and  estimation  of 
background  signal  levels  for  target  detection  sensors.  Direct  modeling  of 
this  radiation  requires  detailed  consideration  of  the  altitude  distribution  of 
infrared  active  molecular  species,  the  line  shape  and  structure  within 
absorption  bands,  and,  in  the  upper  atmosphere,  the  population  and  loss 
mechanisms  of  molecules  in  various  vibrational-rotational  levels.  A 
familiar  operational  computer  code  for  performing  such  direct  radiation 
calculations  is  that  due  to  Degges  (1974).  For  many  applications,  the  accuracy 
of  predictions  based  on  models  is  not  sufficient.  Actual  measured  radiation 
levels  are  required.  On  the  other  hand,  such  measurements  are  expensive, 
difficult  to  obtain,  and  often  not  available  for  the  desired  observation  geometry. 
Under  certain  circumstances  it  is  possible  to  scale  radiation  measurements 
made  for  one  observation  geometry  to  another.  Such  a  scaling  transformation 
is  considered  in  this  paper.  In  particular,  the  problem  addressed  is  that  of 
determining  the  downwelling  atmospheric  radiation  incident  on  an  upsounding 
sensor  from  the  tangent-height  profile  of  earth-limb  radiation.  The  derivation 
of  the  transformation  is  considered  in  Section  2.  An  important  limitation  of 
the  method  is  that  it  can  be  applied  only  under  conditions  where  the  atmosphere 
is  optically  thin.  Consideration  of  this  limitation  is  made  in  Section  3. 

Example  applications  of  the  transformation  method  are  made  for  the  9.  6-pm 
03  and  15-pm  CO^  bands  in  Section  4. 
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2.  Theoretical  Formulation 

a.  General  transformation 

The  fundamental  function  describing  the  radiative  characteristics 
of  a  gaseous  medium  such  as  the  atmosphere  is  the  volume  emission 
coefficient  J.  This  function  describes  the  energy  (e.g.  ,  photons  or  joules) 
emitted  per  unit  time  per  unit  volume  of  atmosphere.  In  general.  J  is  a 
function  of  position  in  the  atmosphere  and  direction  of  emission.  For  the 
present  application,  it  will  be  assumed  that  J  is  a  function  only  of  altitude 
above  the  earth's  surface  (and  not  of  geographic  position  on  the  earth's  surface) 
and  is  isotropic.  Also,  J  is,  in  general,  a  spectroscopic  parameter. 

In  the  present  context,  it  is  assumed  that  J  is  the  mean  value  in  some  sensor 
bandpass.  In  the  development  presented  here,  the  spectral  resolution  of  the 
sensor  is  immaterial  and  may  range  from  essentially  monochromatic  to 
extremely  wide  band.  Radiation  emitted  along  some  line  of  sight  through  the 
atmosphere  is  determined  as  an  integral  of  J  along  the  line  of  sight.  The 
geometries  relevant  to  the  present  problem  are  shown  in  Figs,  la  and  2a. 

(Note  that  both  altitude  and  tangent  height  are  measured  from  earth  center. 

R  is  the  earth  radius. )  Figure  la  illustrates  the  tangent  viewing  geometry 
for  limb  scanning;  z  is  the  tangent  height  and  N^(z)  is  the  in-band  radiance 
emitted  from  the  line  of  sight  s  at  z.  Figure  lb  is  a  qualitative  representa¬ 
tion  of  the  profile  of  N^(z);  H  and  are  the  upper  and  lower  limb  scan 
bounds.  Figure  2a  illustrates  the  upsounding  viewing  geometry,  h  is  the 
sensor  altitude  and  0  is  the  local  zenith  angle  of  observation;  N^h,  0)  is  the 
radiance  presented  to  the  sensor.  Figure  2b  is  a  qualitative  representation 
of  the  profile  Ny(h)  for  some  fixed  0.  The  limb  scanning  and  upsounding 
sensors  are  assumed  to  have  the  same  bandpass  and  fields  of  view.  In  terms 
of  the  volume  emission  function,  the  radiative  transfer  equations  for  these 
two  lines  of  sight  are 


l 

r 
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and 


H 


Nl(z)  =  2  J  J(r) 


rdr 


H 


N^h,  0)  =  [  J(r) 
h 


(  r2  -  z2) 


rdr 


(  r2  -  h2sin2e) 


(1) 


(2) 


These  equations  are  well-founded  in  the  general  theory  of  radiative  transfer 
in  cylindrical  geometry,  but  are  based  on  three,  as  yet  undiscussed, 
assumptions.  The  first  is  that  the  fields  of  view  of  the  two  sensors  are  infini¬ 
tesimally  small.  The  second  is  that  H  is  large  enough  for  us  to  approxi¬ 
mate  J(r)  =  0  for  r  s  H.  The  third,  and  most  important,  is  that  the  radiation 
processes  are  occurring  in  the  optically  thin  regime.  This  last  assumption 
is  treated  in  Section  3. 

Equations  (1)  and  (2)  form  the  basis  of  the  development.  The  goal 
is  to  relate  Ny(h,  6)  directly  to  NL(z).  The  method  is  the  elimination  of 
J(r)  between  (1)  and  (2).  Equation  (1)  is  a  form  of  Abel's  integral 
equation  in  which  NT  (z)  is  considered  as  a  known  function  and  J(r)  is  the 
unknown  function.  The  solution  for  J(r)  is  given  by  the  Abel  transformation 


J(r) 


H 


r 


dNL(z) 

dz 


dz 


I 


Substitution  of  this  result  into  (2)  yields 


H  H 

Nyfh.G)  =  Jj  ^  F(r,  z)dzdr 
r=h  z=r 


(3) 


(4) 
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where 


F(r,  z)  =  -  f 


i  dNL<*> 


TT 


dz 


(*2-2f(-2 


z  -  r  )  r  -  h  sin  0 , 


T 


(5) 


1  quation  (4)  is  a  simple  area  integral,  and  the  order  of  integration  over 
r  and  z  can  be  interchanged  to  give  the  equivalent  integral 


Ny(h,  0)  = 


H 

ft 

z 

N  (h,  0)  =  l 

j 

z=h 

r=h 

)  into  (6)  yields 

H 

z 

1  (  dNL<*> 

f  . 

ir  J  dz  | 

J  , 

F(r,  z) drdz. 


(6) 


rdr 


T  /  2  2\*  (  2  Z  .  2  \ 

h  |z  -  r  J  |r  -  h  sin  0^ 


f 


dz.  (7) 


The  inner  integral  is  available  in  tables  of  standard  integrals  and  gives 
the  result 


H 

Vh-  e>-7  J 


dNL<*»  „  -./V-h2 

dz  ten 


(8) 


Finally,  an  integration  by  parts  with  u  =  tan"  [(z2  -  h2)/(h2  -  d2)] 
and  dv  =  [dN^(z)/dz]  dz  yields  the  final  transformation  as 


H 

Ny(h,  0)  =  J  w(h,  0,  z)NL(z)dz 


(9) 


where 


w(h,  0,  z)  = 


hcos  0 


(2  .2  .  2  _ \  /  2  Til 

(  z  -  h  sin  0 1  I  z  -  h  ) 


(10) 
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These  two  equations  comprise  the  desired  transformation.  The  upsounding 
radiance  profile  Nu(h,  0)  is  determined  directly  as  a  weighted  integral 
over  the  limb  scanning  profile  NL(z).  Example  applications  of  this  trans¬ 
formation  are  given  in  Section  4. 

The  most  significant  feature  of  the  transformation  is  that  no  reference 
is  made  to  the  emission  function  J(r).  If  it  were  necessary  for  J(r)  to  be 
known,  it  would  either  have  to  be  calculated  or  measured.  The  difficulty  of 
calculating  J(r)  was  discussed  briefly  in  the  introduction.  Measurement 
of  J  can  be  made  by  inversion  analysis  (Gille  and  House,  1971;  Russell  and 
Drayson,  1972).  For  the  optically  thin  case,  (3)  is  the  inversion.  An 
unavoidable  aspect  of  such  inversion  analysis  is  the  propagation  of  experi¬ 
mental  noise.  For  example,  in  (3),  even  the  slightest  of  random  fluctuations 
(noise  fluctuations,  not  structural  variations)  in  the  measured  profile  of 
N  (z)  can  propagate  through  the  inversion  to  yield  extremely  noisy  results 
for  J(r).  Great  care  and  some  form  of  data  smoothing  is  generally  required 
to  perform  meaningful  inversions.  The  principal  source  of  noise  amplifi¬ 
cation  in  the  inversion  of  (3)  is  the  need  to  consider  the  derivative  of  N^(z). 
In  the  transformation  procedure  considered  here,  this  derivative  is  eliminated 
in  the  final  integration  by  parts  leading  to  (9).  [Note:  an  integration  by  parts 
cannot  be  used  to  eliminate  the  derivative  in  (3)  ]. 

The  weighting  function  w(h,  0,  z)  is  shown  in  Fig.  3  for  two  zenith 
angles.  The  important  feature  to  note  is  that  the  upsounding  radiance  at 
altitude  h  is  determined  primarily  by  the  magnitude  of  the  limb  profile  in 
the  vicinity  of  z  =  h.  In  fact,  as  0-  it/ 2,  w(h,  0,  z)  approaches  a  delta 
function  and  yields  the  obvious  result  N^(h,  tt/2)  =  ^  N^(h)  • 
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where 


and 


where 


F.  =  f..  .  -  f. 
1  l+l  1 


*i  -  hcosG 


1  -1 
tan 


(«f  -  ^ 


hcosG 


(14) 


Gi  '  8i+l  "  gi 


z  /z2  -  h2'* 

+  j)  tan0 

.  -1 
sin 

h  +  z.sinG 

l 

.  -1 
-  sin 

h  -  z.sinG 

l 

zi  lzi  h  / 

z.  +  hsinG 

z.  -  hsinG 

. 

l 

l 

. 

Substitution  of  (14)  and  (15)  into  (13)  and  a  collection  of  terms  on 
Nt  (z.)  yields  the  final  quadrature  formula 

Lj  1 


N+l 


Nu«>  e>  *  ^  X)  VjNL(i) 


1=J 


where 


^  -6i,N+l> 


+  (1  -  6.  .) 


'  i+l  VV1  +Aj/  ST 

(gi-gi.l) 


(f .  -  f.  . ) 
i  i-l' 


A. 

l 


Zi- 1  A. 


(16) 


and  where  6„„,  =  0  if  n  /  m  and  =  1  if  n  =  m. 
nm 


(15) 
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c. 


Exponential  a 


* }  # ' 


jjL 


St 


In  an  isothermal,  exponential  pressure  atmosphere  with  a  constant 
species  mixing  ratio,  the  source  function  is  of  the  form  J(r)  “•  Jq  exp(-r/cr) 
where  a  is  the  exponential  pressure  scale  height.  Use  of  this  form  in 
(1)  gives  (with  H  -*  oo  ) 


Nl(z)  =  2  JQ  z  Kjfz/er)  (17) 

where  Kj  is  the  modified  Bessel  function  of  order  1.  The  argument 
z/ar  of  Kj  is  large  (since  ziR*1  6371km  and  a  ~7km),  and  the 
asymptotic  limit  of  can  be  used  to  obtain 

NL(z)  ~  J0(2ircrz)^  e“Z^a.  (18) 

The  tangent  height  variation  predicted  by  (18)  is  very  nearly  exponential 
(with  the  same  scale  height  as  J)  since  the  preexponential  variation  with 
z  is  very  small  within  the  atmosphere  (i.  e. ,  0  sz'  =  z  -  R  ^  500  km).  A 
very  nearly  exponential  fall  of  NJz)  above  z  -  R  150km  is  a  common 
prediction  of  the  Degges  model. 

If  an  exponential  profile  is  assumed  for  N^(z),  then  a  very  simple 
approximate  transformation  can  be  derived  from  the  general  transformation 
of  (9).  Let 

NL(«)  =N0e-i/"r' 


14 


(19) 


with  H  -♦  oo  .  Substitution  of  this  profile  into  (9)  and  the  change  of 
variable 


X  = 


(20) 


yields 


Nu(h,  0)  =  Nq  cos  0 

TT 


J 


expf -(h/ort)  (1  +X2)»1 
X2  +  cos^  0 


dX  . 


(21) 


For  typical  a',  h/or'  >>  1  and  the  integrand  of  (21)  is  very  sharply 
peaked  about  X  =  0.  Most  of  the  contribution  to  the  integral  will  come 
from  regions  where  X  is  small.  To  a  high  degree  of  accuracy,  then, 
the  argument  function  \1  +  *  )  may  be  replaced  with  1  +  X  /2.  Doing  so, 
and  evaluating  the  resulting  standard  integral  yields 


Nu(h,  0)  -  *NL(h)  f(P) 


(22) 


where 

2 

f(P)  =  eP  [1  -  erf(P)], 


P  = 


cos  8 
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and  where  (19)  has  been  used  to  express  the  final  result  in  terms  of 
N^(h).  Note  that  the  approximation  has  combined  the  two  independent 
parameters  h  and  8  into  the  single  parameter  (3  .  In  the  limit  (3  -»  0 
(e.g. ,  9  -*  ir/2),  f({3)  "*  1  and  NTT(h)  -*  £  NT  (h),  which  is  the  correct  result 
In  the  limit  (3 -*  oo  (e.g. ,  h/ed  large  and  8  -•  0) ,  f(P)**  l/P**  and 
N^h)-*  N^(h)/[(  2  irh/or')^cos  8]  . 

This  approximation  is  useful  even  if  the  limb  profile  is  not  exactly 

exponential,  particularly  for  ©  near  tt/2.  A  qualitative  measure  of  the 

extent  of  tangent  height  over  which  NT  would  have  to  be  reasonably 

exponential  in  order  for  the  approximation  to  apply  is  the  half  width  of  the 

2  2 

weighting  function  cos0/(X  +  cos  8).  This  width  is  X^  =  cos  0.  The 
corresponding  width  in  z  about  h  is 

Az  -  h(l  +  Xt2)*-  h 


Az  is  tabulated  in  Table  1  for  0  £  0  £  ir/2.  As  8  approaches  tt/2,  the 
altitude  region  over  which  Nj^  would  have  to  be  exponential  in  order  for 
the  approximation  to  apply  becomes  small.  The  accuracy  of  the  approxi¬ 
mation  for  the  worst  case  of  application,  that  is,  8  =  0,  is  illustrated  in 
Section  4. 


Table  1.  Width  Around  h  Over  Which  Should  Be 
Approximately  Exponential  for  Eq.  (22) 
To  Apply. 


0  (degrees)  Jz(km) 


0 

2680 

30 

2089 

60 

763 

70 

368 

80 

97 

85 

25 

87 

9 

89 

1 

90 

0 

17 


3.  Optical  Thinness 


The  transformation  of  limb  data  to  upsounding  data  requires  the 
assumption  that  the  atmosphere  is  optically  thin  in  the  sensor  bandpass 
for  the  lowest  tangent  height  scanned.  The  most  extreme  criterion  for 
thinness  is  that  the  absorptance  over  the  tangent  line  of  sight  at  the  center 
of  the  strongest  spectral  line  in  the  bandpass  be  small  (or^.  <  0.  1  or  0.  9). 

For  a  tangent  height  z,  the  center  line  absorptance  is 


«c=  1  -  exP 


oo 

i 


-2  |  c(r)  p(r)  kQ(r)  rdr 


(*2-  4 


(24) 


where  c  is  mole  fraction  of  active  species,  p  is  total  atmospheric  pressure,  and 
kQ  is  the  absorption  coefficient  at  the  line  center,  which  for  a  general  Voigt 
line  is 


where 


k<r)=  /inA*  8M  K0[a(r)] 

V  "  /  YD<r> 

KQ(a)  =  ea  ^1  -  erf  (a)j  , 


(25) 


(26) 
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n 


a(r)  = 


(Sn  2) 


Yt  (r) 


i  yl 


Yn<r) 


(27) 


and 


YL(r)  =  Y0P(r)  . 


(28) 


S(r),  Yjj(*')  and  YD(r)  are,  respectively,  the  line  strength,  the  pressure 
broadening  (Lorentz)  component  of  the  line  width,  and  the  Doppler  broaden¬ 
ing  component  of  the  line  width.  Near  the  centers  of  the  strong  molecular 
absorption  bands,  the  line  strengths  are  reasonably  independent  of  tempera¬ 
ture  (and  thus  of  altitude)  and  we  take  S(r)  **S.  Also  assume  c(r)  •*  c  and 
Yp(r)  Yd  .  Then,  for  an  exponential  pressure  atmosphere 


P(r)  =  p0e"(r'R)/a  (29) 

with  scale  height  or  ,  the  absorptance  at  tangent  height  z  can  be  expressed 
by 


<*C(P)  =  1  -  exp 


(30) 


where  (3  =  z/a  and 
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F(p)  =  P  j 


KM 


a0  ' 


w 


a0  =  . 


p  =  r/z,  and  Pn  =  R/or.  F  was  evaluated  numerically  with  R  =  6371  km, 

u  _  j  _i 

or  =  6.65  km,  =  1  atm,  y^  =  0.08  cm  /atm,  and  y^  =  0.0005  cm 
(thus  »q  =  133  and  p^  =  958).  The  result  is  shown  in  Fig.  4.  Also  shown 
is  the  function  Fp(z)  that  would  obtain  for  a  pure  Doppler  line.  This  latter 
case  is  obtained  from  (31)  in  the  limit  a_  —  0  and  is 


oo 

FD(P)=  p| 


-<pp-p0>  pd 


fpM5 


Kj  (P) 


where  Kj  is  the  modified  Bessel  function  of  order  1.  Since  P  is  large 
(P  a  Pq  =  958),  asymptotic  limit  Kj  (x  )  ~  (n72x)^  e  *  can  be  used  to 
obtain 

1  -(P-P0) 

Fd(P)=  (,p/Z)*e  V  (3 


The  comparison  in  Fig.  4  confirms  the  well-known  result  that,  for  all 
practical  purposes,  a  purely  Doppler  line  shape  can  be  assumed  for  altitudes 
above  ~  50  km. 

Equation  (30)  holds  for  an  isolated  line.  The  most  dense  crowding  of 
spectral  lines  occurs  in  the  Q  branch  of  the  15-p.m  CC>2  band.  To  a  good 


approximation,  the  line  spacing  in  the  Q  branch  is  given  by 


6  -  (2  J  +  1)  (b'  -  b") 


(34) 


where  J  is  the  rotational  quantum  number  and  B  and  b'  are  the  rotational 
constants  of  the  upper  and  lower  vibrational  levels,  respectively.  The  most 
dense  crowding  occurs  at  J  =  0,  and  with  b'  -  b"  =  10.  2  x  10"4cm"1 
(Kaplan  and  Eggers,  1956),  (34)  yields 


6  . 
min 


0.  001  cm" 


The  Doppler  line  width  is  yD  0.  0005  cm-1  and  is  half  this  worst  estimate 
of  line  spacing.  The  Lorentz  half  width  is 


VL  “*  Y0  P0 


e-(z-R)/c^ 


With  Vq  =  0.08  cm  *,  p^  =  1  atm,and  or  =  6.65km,  the  Lorentz  width  will  be 

smaller  than  6  .  for  altitudes  z-R>  29km.  Thus,  the  assumption  of 
mm  r 

isolated  lines  even  in  the  Q  branch  is  justified  for  altitudes  above  ~  30  km. 

The  line  center  absorptance  or  was  computed  for  each  of  the  major 
absorbing  atmospheric  species.  These  species  and  the  locations  of  their 
absorption  bands  are  tabulated  in  Table  2.  Within  each  band,  the  strongest 
line  was  located  by  searching  the  AFGL  line  atlas  (McClatchey  etal.  ,  1973). 
This  strength  and  the  location  v  of  the  line  are  also  listed  in  Table  II.  An 
estimate  of  the  likely  maximum  value  of  species  concentration  near  100km 
altitude  is  also  given  in  the  table  (Degges,  1974).  From  these  data,  the 
minimum  altitudes  for  which  a  _  £  0.  1  were  computed  and  are  listed  in  the 

^  I 

final  column  of  the  table.  A  plot  of  =  1  -  a  ^  versus  z  =  z-R  for 

CC>2  is  shown  in  Fig.  5.  Except  for  the  CO^  band,  the  minimum  altitudes 

are  ~100  _+  10km.  For  CC>2»  the  altitude  is  somewhat  higher,  at  about  135km. 
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Table  2.  Maximum  Line  Strength  and  Concentration  for 
Important  Atmospheric  Absorption  Band 
(5  to  20  pm) 


SpVC  It*  Si 

ILmd  (jimi 

S  (cm‘^/atm  at  296  K) 

max 

(cm*  1  ( 

dmolf  fraction  1 
at  ar  J  00  km 

7.  .  (km) 

min' 

II  ,o 

7.  94 

168S 

3  (-61 

108 

ch4 

7.  7 

1.83 

1  333 

1  (-61 

91 

n^o 

7.  8 

4.  36 

1  £98 

1  (-6) 

J01 

».  6 

0.  99 

I  OS  3 

l  1-61 

91 

CO  , 

1  * 

7.  69 

b67.  7 

i  (—43 

136 
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Because  of  the  number  of  assumptions  made  in  the  derivation  of  (30) 
and  because  of  the  quality  of  data  (particularly  c  and  or )  used,  these  altitude 
results  should  be  viewed  only  as  rough  upper  estimates  of  the  minimum 
altitude  at  which  optical  thinness  obtains.  For  example,  results  obtained 
with  the  Deggqs  code  for  the  13-  to  16-p.m  CO^  bandpass  (see  next  section) 
indicate  that  optical  thinness  prevails  at  lower  altitudes.  Under  optically  thin 
conditions,  the  "upsounding"  radiance  at  9  =  ir/2  should  be  exactly  half  the 
limb  radiance,  and  the  function 


2Nu(z,  rr/2) 

N^)~ 


(35) 


should  be  unity.  The  function  actually  deduced  from  the  code  calculations  is 

I 

shown  in  Fig.  5  and  displays  optical  thinness  (i.  e. ,  p  >0.9)  down  to  z  =115  km. 


In  summary,  near  the  centers 
mation  method  should  not  be  applied 


of  strong  absorption  bands  the  transfor- 

I 

for  tangent  heights  lower  than  z  ~  100km. 
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4.  Example  Application 


Synthetic  limb  radiance  profiles  for  z  2  100  km  were  generated  with 
the  Deggs  code  for  CC>2  in  the  13-  to  16-|i.m  bandpass  and  for  in  the 
9-  to  10-pm  bandpass.  The  CO^  profile  is  shown  in  Fig.  6  for  both 
nighttime  and  daytime  conditions.  Upsounding  profiles  for  9  =  0,  60,  80 
and  90°  were  also  computed  with  the  code  for  the  nighttime  condition  and 
are  shown  in  Fig.  7.  Results  of  the  application  of  the  numerical  transforma¬ 
tion  and  the  exponential  approximation  to  the  night  limb  profile  are  shown  in 
Fig.  8  for  0  =  0°.  Except  in  the  region  from  100  to  ^  130km,  the  numerical 
transform  reproduces  the  model  prediction  exactly.  The  discrepancy  below 
130km  is  a  result  of  the  violation  of  the  assumption  of  optical  thinness  explicit 
in  the  formulation  of  the  transformation.  For  other  zenith  angles,  the  agree¬ 
ment  between  the  code  and  numerical  transform  predictions  is  similar.  The 
agreement  for  9  =  0°  between  the  code  prediction  and  the  exponential  profile 
approximation  is  essentially  exact  for  altitudes  above  300km.  This  is  to  be 
expected  since  above  this  tangent  height,  the  fall-off  of  the  limb  profile 
(Fig.  6)  is  almost  exactly  exponential.  The  error  of  the  approximation  is 
greatest  (not  counting  the  error  due  to  the  violation  of  the  assumption  of 
optical  thinness)  at  130km  and  amounts  to  ~  16%.  The  case  0  =  0°  is  the 
worst  condition  of  application  of  the  exponential  approximation,  and  the  maxi¬ 
mum  error  decreases  monotonically  to  zero  as  9  goes  to  90°. 

Comparisons  between  code  predictions  and  transformation  results  for 
are  shown  in  Figs.  9-11.  The  atmosphere  is  optically  thin  down  to  about 
90km  in  this  spectral  region,  and  the  agreement  between  the  code  prediction 
and  the  numerical  transformation  result  is  essentially  exact  for  all  altitudes 
and  zenith  angles  (see  Fig.  11  for  0  =  0°).  Below  and  above  the  sharp  transition 
in  the  limb  profile  (Fig.  9)  at  120  km,  the  exponential  approximation  is  also 
quite  accurate.  In  the  transition  region,  however,  the  exponential  profile 
approximation  underestimates  the  upsounding  profile  by  as  much  as  a  factor 
of  two. 
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Fig.  7.  Code  Predictions  for  CC>2  Upsounding  Profiles 
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Fig.  8.  Comparisons  Among  CO.  Upsounding  Profile 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  Investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  nllltary  concepts  and  aystena.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laborato¬ 
ry  personnel  In  dealing  with  the  many  problems  encountered  In  the  Nation's 
rapidly  developing  apace  aystasM.  Expertise  In  the  latest  scientific  develop¬ 
ments  la  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aetophyelce  Laboratory:  Aerodynamics;  fluid  dynamics;  plasmadynamlcs; 
chemical  kinetics;  engineering  mechanics;  flight  dynamics;  heat  transfer; 
high- power  gas  lasers,  continuous  and  pulsed,  IR,  visible,  UV;  laser  physics; 
laser  resonator  optics;  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  optical  hack- 
g rounds;  radiative  transfer  and  atmospheric  transmission;  thernal  and  state- 
apaclflc  reaction  rates  In  rocket  plumes;  chemical  thermodynamics  and  propul¬ 
sion  chemistry;  laser  Isotope  separatism;  chemistry  and  physics  of  particles; 
apace  environmental  and  contamlnat Ion  affects  on  spacecraft  oaterlals;  lubrica¬ 
tion;  surface  chemistry  of  Insulators  and  conductors;  cathode  na ter la Is;  sen¬ 
sor  materials  and  sensor  optics;  applied  laser  spectroscopy;  atonic  frequency 
standards i  pollution  and  toxic  materials  monitoring. 

Electronics  Research  Laboratory:  Electromagnetic  theory  and  propagation 
phenomena;  microwave  and  semlcnndertor  devices  and  Integrated  circuits;  quan¬ 
tum  electronics,  lasers,  and  electro-optica;  conmunlcatlon  sciences,  applied 
electronics,  superconducting  and  electronic  device  physics;  millimeter— wave 
and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  nan  oaterlals;  composite 
materials;  graphite  and  ceramics;  polymeric  materials;  weapons  effects  and 
hardened  materials;  materials  for  electronic  devices;  dimensionally  stable 
materials;  chemical  and  structural  analyses;  stress  corrosion;  fatigue  of 

metals. 

Space  Sciences  laboratory:  Atmospheric  and  Ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  tho  atmosphere,  aurorae 
and  alrglow;  nsgnetospherlc  physics,  eoomlc  rays,  gsnaratleo  and  propagation 
of  plasma  moves  la  tho  aagnetosphsro;  solar  physics,  x-ray  astronomy;  the  offsets 
of  nuclear  explosions,  magnetic  stoma,  and  solar  activity  on  tho  north's 
ataosphere.  Ionosphere,  sad  aagnatoaphare;  tho  offsets  of  optical,  electromag¬ 
netic,  and  particulate  radiations  In  space  on  space  systems. 


